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ON THE NEUTRINO FLUX FROM GAMMA RAY BURSTS: AN OVERVIEW 



DAFNE GUETTA 
Osservatorio astrofisico di Arcetri, Largo E. Fermi 5, 
50100 Firenze, Italy 



Observations imply that gamma-ray bursts (GRBs) are produced by the dissipation of the 
kinetic energy of a highly relativistic fireball. Photo-meson interactions of protons with 7- 
rays within the fireball dissipation region are expected to convert a significant fraction of the 
fireball energy into high energy neutrinos. In this talk we summarize the results of an analysis 
of the internal shock model of GRBs, where production of synchrotron photons and photo- 
meson neutrinos with energies ~ 10 14 eV are self-consistently calculated. These neutrinos will 
be coincident with the GRB. We show that the fraction of fireball energy converted into high 
energy neutrinos is not sensitive to uncertainties in the fireball model parameters, such as 
the expansion Lorentz factor and characteristic variability time. Other processes of neutrinos 
emission from GRBs are also reviewed. Photomeson interactions within the plasma region 
shocked by the reverse shock, may produce a burst of ~ 10 18 eV neutrinos following the GRB 
on the time scale of 10 s. Inelastic p-n nuclear collisions result in the production of a burst 
of ~ 10 GeV neutrinos in coincidence with the GRB. Planned 1 km 3 neutrino telescopes are 
expected to detect ten 100 TeV neutrino events and several 10 18 eV events, correlated with 
GRBs, per year. A suitably densely spaced detector may allow the detection of several 10 
GeV events per year. 



1 Introduction 

The characteristics of 7-ray bursts (GRBs), bursts of 0.1 MeV-1 MeV photons lasting for a 
few seconds §, suggest that the observed radiation is produced by the dissipation of the kinetic 
energy of a relativistically expanding wind, a "fireball," at cosmological distance (see, e.g., El for 
review) . The recent detection of delayed low energy (X-ray to radio) emission (afterglow) from 
GRB sources (seeH for review), confirmed both the cosmological origin of the bursts and the 
standard model predictions for afterglows, that result from the collision of an expanding fireball 
with its surrounding medium (seeu for review). 

Within the fireball model framework, observed 7-rays are produced by synchrotron emission 



of electrons accelerated in the fireball dissipation region. In this region the plasma parameters 
allow Fermi shock acceleration of protons up to energies > 10 20 evi'l; (seei for a recent review). 

High energy protons may photo-interact with the fireball photons leading to emission of 
high energy neutrinos!!!. Electrons and protons can be accelerated to high energy by internal 
shocks within the expanding wind or by the reverse shock during the transition of the plasma to 
self-similarity. In the first case a burst of ~ 100 TeV neutrinos is expected in coincidence with 
the GRB, and if protons are accelerated in the reverse shock a burst of ~ 10 18 eV neutrinos is 
expected to follow the GRB pji r a 1 time scale of 10 s. Lower energy neutrinos may be produced 
by inelastic nuclear collisions l 10 H n l. 

In this review new results found in Guetta, Spada & Waxman 2001 (GSWOl)ll! on the v 
flux expected in the internal shock model scenario are presented and recent work on the high 
energy neutrino production by GRB fireballs is summarized. Implications for planned ~ 1km 3 
neutrino telescopes (the ICECUBE extension of AMANDA, ANTARES, NESTOR; seelll for 
review) are discussed. High energy neutrino production in internal shocks is discussed in § g, 
other mechanisms for neutrino production in GRBs are summarized in § || and the implications 
for the planned neutrino telescopes are presented in § |j. 



2 Internal shocks neutrinos 



High energy neutrinos are produced by the decay of charged pions, ir + — > fi + + 
v e + Vu -\-Vp,i created in interactions between the fireball photons and accelerated protonsE 
The characteristic energy of these neutrinos is determined by the relation between the observed 
photon energy, E^, and the accelerated proton's energy, E p , at the photo- meson threshold of 
the A-resonance. In the observer frame, 

E y E p = 0.2GeV 2 r 2 . (1) 

For typical observed 7-ray energy of 1 MeV and Lorentz factors of the expanding fireball T ~ 300, 
proton energies i? p ~2x 10 7 GeV are required to produce neutrinos from pion decay. Typically 
the neutrinos receive 5% of proton energy leading to v of ~ 10 15 eV. 

The predicted flux of ^ 10 14 eV neutrinos produced in internal shocks is determined by 
the fraction /„- of fireball proton energy lost to pion production. This fraction is given by 
f w = min(l, At/t n ) where At is the comoving shell expansion time and t 7r is_the proton photo- 
pion energy loss time. The value of this fraction has been estimated to be&E-^ 

U(E P ) n0.2wm(l, E p /E*) ^> 52 . (2) 

The 7-ray luminosity L 7 , the photon spectral break energy E^ (where luminosity per logarithmic 
photon energy interval peaks), the observed variability time At and the wind Lorenz factor T 
are normalized in Eq. |2]to their typical values inferred from observations: L 7 = 10 52 erg/s, £* = 
\E MeV MeV, At = 10 _2 s, and T = 10 2 ' 5 . The proton break energy, E p , is the threshold proton 
energy for interaction with photons of observed energy E*, E p = (2/E^ Mev)^i.5 x 10 7 GeV. For 
a typical burst, E ~ 10 53 erg at redshift z ~ 1, Eq. |2| implies a neutrino fluence 

where the neutrino break energy is given by El m 10 15 r| 5 /^ MeV eV, whereas a = for 
E v > E h v and a = 1 for E v < E b v . 

The value of f n , Ecl ||, is strongly dependent on T. It has recently been pointed out 
by Halzen and Hooper that if the Lorentz factor V varies significantly between bursts, the 



resulting neutrino flux will be dominated by a few neutrino bright bursts, and may significantly 
exceed the flux given by Eq. [||, derived for typical burst parameters T his may strongly enhance 
the detectability of GRB neutrinos by planned neutrino telescopes 113. 

In this section we present the most important results of the work GSWOllil whose main goal 
has been to determine the allowed range of variation in the fraction of fireball energy converted 
to high energy neutrinos, under the assumption that GRBs are produced by internal dissipative 
shocks in an ultra-relativistic wind. 



2.1 Model dynamics and neutrino production 

In the fireball model of GRB i a compact source, of linear scale Rq ~ 10 6 cm, produces a 
wind characterized by an average luminosity L w ~ 10 52 ergs -1 -10 53 ergs -1 and mass loss rate 
M = L w /rjc 2 . Variability of the source on time scale At, resulting in fluctuations in the wind bulk 
Lorentz factor T on a similar time scale, leads to internal shocks that can reconvert a substantial 
part of the kinetic energy to internal energy. It is assumed that this energy is then radiated as 
7-rays by synchrotron (and inverse-Compton) emission of shock-accelerated electrons. 

In GSW01 we have used an approximate model for the unsteady wind following Spada, 
Panaitescu & Meszaros 2000 @, and Guetta, Spada & Waxman 2000 (GSW00) 0. The wind 
flow is approximated as a set of discrete shells. Each shell is characterized by four parameters: 
the ejection time tj, where the subscript j denotes the j-ih shell, the Lorentz factor Tj, the mass 
Mj, and the width A-,-. Since the wind duration, t w ~ 10 s, is much larger than the dynamical 
time of the source, td ~ Rq/c, variability of the wind on a wide range of time scales, td < t v < t w , 
is possible. For simplicity, we considered a case where the wind variability is characterized by 
a single time scale t v , in addition to the dynamical time scale of the source td and to the wind 
duration t w . Thus, we considered shells of initial thickness Aj = ctd = Rq, ejected from the 
source at an average rate t~ . In GSW00 we examined the dependence of the observed 7- 
ray flux and spectrum on wind model parameters, taking into account both synchrotron and 
inverse-Compton emission, and the effect of e^ pair production. We have shown that in order 
to obtain 7-ray fluxes and spectra consistent with observations, large variance is required in the 
wind Lorentz factor distribution. We therefore restricted our study to the bimodal case, where 
Lorentz factors are drawn from a bimodal distribution, J?a — T m or Tj = Tm ~ V* ^ T m , with 
equal probability. The time intervals tj + \ — tj were drawn randomly from a uniform distribution 
with an average value of t v . In GSW00 we considered two qualitatively different scenarios 
for shell masses distribution, shells of either equal mass or equal energy, and concluded that 
observations favor shells of equal mass. We therefore restricted our analysis to the equal shell 
mass case. 

Once the shell parameters had been determined, we calculated the radii where collisions occur 
and determined the photon and neutrino emission from each collision. In each collision a forward 
and a reverse shock are formed, propagating into forward and backward shells respectively. We 
assumed that a fraction e e of the protons shock thermal energy is converted to electrons and a 
fractio €3 to magnetic field. We adopted these fractions to be close to equipartition, e e = 0.45 
and €b ^ 0.01. We further assumed that both electrons and protons are accelerated by the 
shocks to a a power-law distribution, dn a /d , y a tx 7^ p , with p ~ 2 for particle Lorentz factors 
7a,min < la < 7a,max- The maximum Lorenz factor is determined by equating the acceleration 
time, estimated as the Larmor radius divided by c, to the minimum between the dynamical and 
synchrotron cooling time. 

The calculation of the emitted v spectrum and flux was carried out according to the method 
described in GSW01, following Waxman & Bahcall 1997 0. Photo-meson production of low 
energy protons, well below requires interaction with high energy photons, well above Ey 
Such photons may be depleted by pair production. For each collision we found the photon 



energy E^, for which the pair production optical thickness, r 77 , is unity. A large fraction of 
photons of energy exceeding max(£^,m e c 2 ) (measured in the shell frame) will be converted 
to pairs, and hence will not be available for photo-meson interaction, leading to a suppression 
of the neutrino flux at low energies. In order to take this effect into account, we used f n = 
min[l, (At/t n ) min(l, t~J)] for protons interacting with photons which are strongly suppressed 
by pair production. 

Neutrino production is suppressed at high energy, where neutrinos are produced by the 
decay of muons and pionswhpse lifetime t exceeds the characteristic time for energy loss due 
to synchrotron emission &EJIL3. We therefore defined an effective f n , f w e s.i as 4/^ times the 
fraction of pions' energy converted to muon neutrinos. 

In each collision a fraction of the kinetic energy of the colliding shell is converted to a flux of 
photons and neutrinos. The energy that is not lost to photons and neutrinos is converted back 
to kinetic energy by adiabatic expansion of the shell. 

2.2 Main results and conclusions 

The main results of GWS01 are presented in Figures 1 and 2. Figure 1 presents the contour 
plot for the effective f n , f n>e ff. taking e# = 0.01. The values shown are averaged over all internal 
collisions, weighted by the energy converted to high energy protons in each collision. Values of 
/tt are shown for several neutrino energy ranges, using the approximation E v = 0.05E p . The 
region in F m —t v plane where E^ > 0.1 MeV is bound in Figures 1 by the dashed line. The 
dash-dotted line outlines the region in which the fraction of wind energy converted to radiation 
exceeds 2% (higher fraction is obtained for larger F m values). Figure 2 shows the contour plots 
of f n , neglecting the effects on neutrino production of synchrotron losses and pair-production. 
/tt approaches unity at low values of F m for low proton (neutrino) energies, and at intermediate 
values of F m for high proton (neutrino) energies. The energy loss of pions and muons due 
to synchrotron emission reduces the fraction of pion energy converted to neutrinos, and hence 
suppresses the effective value of /„-, to the range of values shown in Figure 1. The suppression 
seen in Figure 2 of f n at low F m for high energy protons is due to the fact that for low values of 
F m most collisions occur at small radii, where synchrotron losses of high energy protons prevent 
acceleration of protons to ultra-high energy 0. Ultra- high energy protons are produced in this 
case only by secondary collisions at large radii, where the photon energy density is already low, 
resulting in smaller values of f n . 

As we can see from Figure 1, for wind parameters consistent with observed GRB charac- 
teristics, the effective value of f w at high neutrino energies, E v > E^, is restricted to values in 
the range of « 10% to sa 30%. The weak dependence of f^ >e s, on wind the model parameters, 
in contrast with the strong dependence implied by Eq. is due to two reasons. First, for low 
values of F and At, where large values of are implied by Eq. ^, only a small fraction of the 
pions' energy is converted to neutrinos at high proton energy due to pion and muon synchrotron 
losses (compare figures 1 and 2). Secondth, the observational constraints imposed by 7-ray 
observations imply that the wind model parameters (F, L, At) are correlated. From Figure 
1 we note, also, that the value of f W:e s. does not significantly exceed 20% also in wind model 
parameter regions which are outside the parameter region implied by observations. Our results 
imply that GRB neutrino flux of individual bursts should correlate mainly with the bursts' 7-ray 
flux. 



3 Neutrinos from other processes in GRBs 
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Figure 1: Contour plots of /^cff., the effective value of (defined as 4f n times the fraction of pions' energy 
converted to muon neutrinos), as function of wind variability time t v and minimum Lorentz factor T m , for 
wind luminosity L w = 10 53 erg/s and ts = 0.01. The four panels correspond to four observed neutrino energy 
bins, clockwise from top left: 10 14 eV < E„ < 10 15 eV, 10 15 eV < E v < 10 16 eV, 10 17 eV < E v < 10 18 eV, 
10 18 eV < E v < 10 19 eV. We have used the approximate relation E v = 0.05-Ej, between neutrino and proton 
energy, and assumed a source redshift z — 1.5. The region in the V m -t v plane where E^ > 0.1 MeV is bound 
by the dashed lines. The dash-dotted lines outline the region in which the fraction of wind energy converted to 
radiation exceeds 2% (a higher fraction is obtained at larger Y m values). 
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Figure 2: Contour plots of f n , obtained when the effects on neutrino production of synchrotron losses and pair- 
production are neglected, for the case shown in Figure 1. 



3. 1 Neutrinos from the reverse shock 



Much higher energy, ;> 10 18 eV, neutrinos may be produced at a later stage, at the onset 
of the fireball interaction with its surrounding medium. During the transition to self-similar 
expansion, which takes place on time scale ~ t w , protons and electrons are accelerated in the 
reverse shocks. The combination of low energy photons, 10 eV-1 keV, and high energy protons 
produces neutrinos of ultra-high energy 10 17 — 10 19 eV via photo-meson interactions, as indicated 
by Eq. [1]!. However, the fraction of energy converted to pions depends upon parameters of the 
surrounding medium. If the density of the surrounding gas is that typical of the interstellar 
medium uq ~ 1cm -3 , the photon emission peaks in the X-ray band, E^^y and the luminosity 
is Lx ~ 10 50 erg/s. Using these parameters in Eq. [j], replacing At with t w and considering that 
the self-similar Lorentz factor of the expanding plasma under these conditions is Tbm ~ 250, we 
nd f n = 0.01. Thus the expected neutrino fluence for a typical burst, E w 10 53 erg at z ~ 1, is 

where E h v ~ 10 17 eV, a = 1/2 for E v < E b v and a = 1 for E v > E h v . If the fireball expands into 
a pre-existing wind, the transition to self-similar behaviour takes place at a radius where the 
wind density is n ~ 10 4 cm -3 . Since Tbm decreases with increasing density, according to Eq. |l| 
protons of energy E p > 10 18 eV lose all their energy to pion production in the wind case, and a 
typical GRB at z ~ 1 is expected to produce a neutrino fluence E 



where again £* ~ 10 17 eV, a = for E v < E b u and a = 1 for E v > E b v . 



3.2 Inelastic p-n collisions 

We briefly discuss how lower energy neutrinos may be produced by inelastic nuclear collisions, 
a detailed analysis of these processes is given inEJ'E3. According to most progenitor scenarios 
neutrons can contribute to the barionic component of the fireball as much as protons. During the 
acceleration phase of the fireball protons and neutrons are coupled by nuclear scattering. As the 
fireball expands neutrons and protons may decouple and relativistic relative velocities may arise 
between them, leading to pion production through unelastic collisions. Since the decuopling 
happens when the wind Lorentz factor is ~ 400, neutrons decouple from the accelerating plasma 
prior to saturation, V = rj, only if r] > 400. Charged pions may decay in ~ 10 GeV neutrinos. 
A typical burst, E = 10 53 erg at z ~ 1, with a significant neutron to proton ratio and r\ > 400 
will produce a fluence of ~ 10 -4 cm -2 of ~ 10 GeV neutrinos. 



4 Implications 

The high energy neutrinos predicted in the dissipative wind model of GRBs may be observed by 
detecting the Cerenkov light emitted by high energy muons produced by neutrino interactions 
below a detector on the surface of the Earth (seecJ for a recent review). 

For the internal shocks neutrinos with energies 10 14 — 10 15 eV, the fluence found in Eq. |3| 
implies a detection probability of ~ 0.03 per burst in a km-cube telescope, corresponding to tens 
of events per years correlated in time and direction with GRBs, given the observed GRB rate of 
~ 10 3 yr _1 . The predicted fluence of neutrinos of 10 17 eV produced during the transition of the 
fireball to self-similarity depends strongly upon the conditions of the surrounding gas. In the 
case of a fireball expanding in a typical interstellar medium Eq. || implies a detection probability 



of ~ 10 -4 ' 5 per burst in a km-cube telescope. This probability is higher, ~ 10 -2 ' 5 , in the case 
of fireball expansion into a pre-existing massive star wind. In this case a km 3 telescope should 
detect several muon induced neutrinos per year. 

Inelastic p-n collisions may produce neutrinos of ~ 10 GeV with a fluence of ~ 10 -4 cm ~ 2 
per burst, due to the decoupling of neutrons and protons during the acceleration phase. This 
process is possible only if large neutron to proton ratio and high, > 400, Lorentz factor are 
assumed. Ten events per year are expected in a km 3 detector. Such events may be detectable 
in a suitably densely spaced detector. 

Detection of high energy neutrinos will test the shock acceleration mechanism and the sug- 
gestion that GRBs are the sources of ultra-high energy protons, since > 10 14 eV (> 10 18 eV) 
neutrino production requires protons of energy > 10 16 eV (> 10 19 eV). Detection of ~ 10 GeV 
neutrinos will constrain the fireball neutron fraction, and hence the GRB progenitor. 
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